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RADIOFREQUENCY  RADIATION  AND  LIVING  TISSUE: 
THEORETICAL  STUDIES 


This  report  is  concerned  with  basic  mechanisms  whereby  radiofrequency 
electromagnetic  radiation  can  interact  with  living  tissue.  Our  analysis  of 
basic  mechanisms  suggests  that  radiofrequency  electromagnetic  radiation  may 
be  significant  in  its  effect  on  living  tissue. 

Part  I  of  this  report  outlines  an  elementary  statistical  physics  model 
of  tissue  viewed  as  a  chemical  system.  Part  II  presents  a  nonquantum  pt\ysics 
analysis  of  macromolecul ar  motion  in  an  electromagnetic  field.  The  final 
section.  Part  III,  provides  a  discussion  and  recommendation  for  further 
analysis  and  experimentation. 

PART  I:  STATISTICAL  PHYSICS  MODEL  OF  TISSUE 
THE  MODEL  TISSUE 

The  approach  taken  here  is  to  model  tissue  as  a  mixture  of  two  molecu¬ 
lar  types,  #1  and  #2.  Type  #1  is  thought  of  as  cell  water,  and  type  #2 
consists  of  biomacromolecules.  Type  #1  molecules  can  exist  in  energy  levels 

Uu»  Uj2»  Uj3,  . ,  Ujl ;  while  type  #2  molecules  can  exist  in  energy 

levels  U21,  U22»  U23,  .  U2M*  These  molecules  are  considered  to  be 

only  weakly  interacting  and  each  is  considered  distinguishable  according  to 
its  spatial  position  or  orientation,  so  that  Maxwell -Boltzmann  distribution 
theory  is  valid  (13). 

Let  the  symbol  x^j  represent  the  number  of  molecules  of  type  i  in 
energy  state  U-f j .  If  the  total  number  of  molecules  of  type  #1  is  Nj,  then 

L 

l  *1J  "  "1  (1) 

j-1 


and,  if  the  total  number  of  molecules  of  type  #2  is  N2,  then 
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M 

l  *2j  =  n2  <2) 

J-l 

Each  energy,  level  Ufj  has  an  associated  set  of  gjj  elementary  eigenstates,  as 
schematically  represented  in  Figure  1  (13).  Therefore,  for  the  ensemble  of 
type  #1  molecules,  the  total  number  of  possible  eigenstates  associated  with 
the  occupancy  numbers  xjj,  ....  xjl  is  Gj,  where 


Hi! 

Xu!  — x1Ll 


Likewise,  for  the  type  #2  molecules,  G2  is 


(4) 

(4.1) 


This  value  pfj  can  be  considered  as  the  a  priori  probability  of  the  energy 
level  Ujj.  Substituting  In  the  above  relation,  the  following  equations  are 
obtained: 


Nj! 

n*ll 

PX1L 

*11 

!...x1L! 

vn 

P1L 

M2! 

d*21 

p*2M 

*21 

! . . .  x2m ! 

*21  ' 

P2M 

where  Pj  and  Pg  are  the  probabilities  of  the  occupancy  configurations  x^, 
....  xjl  and  X2i,  ....  X2M»  respectively.  Whence  the  total  system 
probability  can  be  written  as 

P  =  Pi P2 
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Figure  1.  Energy  state  structure  of  the  model  tissue. 


It  Is  shown  in  Appendix  A  that,  if  the  mixture  is  in  a  state  of  equi¬ 
librium,  the  following  relationships  hold  for  njj  (the  most  likely  value  of 
xjj)  and  for  n2j  (the  most  likely  value  of  X2j): 

L  -1 

"li  =  PH  "l  “p  (VkT’  UlPlreX|>,-VkT,l 

"2J  •  V* exp  ‘  VkT) p*j“p,VH,»‘1 

r=l 


In  these  last  two  equations,  k  is  the  Boltzmann  constant  and  T  is  temperature 
given  in  degrees  Kelvin.  Note  that  the  total  system  energy,  U,  can  be  calcu¬ 
lated  using  the  following  equation: 


L 

l 

j=l 


n  U 

lj  lj 


M 

l 

j=l 


n  U 
2j  2j 


=  U 


(6) 


MODEL  OF  RADIATION-INDUCED  TRANSITIONS 

For  purposes  of  illustration,  a  tissue  with  only  four  energy  states  is 
considered.  Type  #1  molecules  have  energy  states  and  Ui2»  and  type  #2 
molecules  have  energy  states  U21  and  U22-  Since  type  #2  molecules  represent 
biomacromolecules,  photon-stimulated  transitions  between  U2i  and  U22  are 
examined. 

Consider  first  the  transitions  between  U21  and  U22  when  the  tissue  is 
not  receiving  externally  applied  radiation  and  is  in  a  state  of  equilibrium. 
If  C21  is  the  equilibrium  rate  of  transition  from  U21  to  U22  per  second  per 
molecule,  and  C22  is  the  analogous  rate  for  the  reverse  transition,  equilib¬ 
rium  requires: 

c21n21  =  c22n22  <7) 

To  analyze  transitions  due  to  radiation,  transition  probabilities  are  used  as 
fo".  ows:  (1)  B21  Is  the  probability  of  transition  from  U21  to  U22  ty  photon 
absorption  per  molecule  in  U21  and  per  unit  field  intensity  (e.g.,  1  mW/cm^); 
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(2)  B22  Is  the  probability  of  transition  from  U22  to  U21  by  stimulated  photon 
emission  per  molecule  in  U22*  and  per  unit  field  intensity. 

Considering  the  case  wherein  the  added  radiation  does  not  significantly 
perturb  mixture  total  energy,  when  the  mixture  achieves  a  new  equilibrium 
with  the  radiation,  the  following  equality  must  hold: 

(B21I  +  C2 1 ) < n2 1  -  a)  *  {B22I  +  C22Hn22  +  4)  (8) 

In  this  equation,  n2i  and  022  are  the  equilibrium  values  calculated  by  equa¬ 
tions  4  and  5,  a  is  the  change  in  021  and  n22  due  to  the  externally  applied 
radiation,  and  I  is  the  radiation  exposure  received  from  an  outside  source  or 
sources.  It  is  well  known  that  B21  =  B22  =  B.  Using  this,  it  is  found  that: 

A  •  (9) 

2BI  +  C21  +  C22 

Equations  4,  5,  6,  and  9  allow  exploration  of  radiation  effects  in  a 
tissue  when  tissue  temperature  does  not  change  significantly.  A  simple 
example  is  calculated  in  the  next  subsection. 

A  CALCULATED  EXAMPLE 

Suppose  the  tissue  is  described  by  parameters  as  given  in  Table  1. 
Recall  that  type  #1  molecules  represent  cell  water  molecules  while  type  #2 
molecules  represent  biomacromolecules.  In  Table  1,  the  value  used  for  Nj  is 
the  number  of  water  molecules  in  a  sphere  with  a  10-pm  radius,  the  sphere 
representing  a  small  cell.  The  value  for  N2  is  an  estimated  number  of  enzyme 
molecules  of  a  particular  type  in  a  single  cell.  The  energy  levels  in  the 
enzyme  are  chosen  so  that  transitions  can  be  caused  by  500-MHz  radiation. 

The  upper  enzyme  mode  with  energy  U22  is  given  a  low  probability,  P22»  to 
achieve  a  mode  that  is  poorly  populated  in  the  absence  of  external  radiation. 
For  example,  the  energy  U22  roay  correspond  to  a  first  step  in  enzyme  dena- 
turation,  which  may  not  be  favored  in  unirradiated  tissue. 


TABLE  1.  TISSUE  SPECIFICATIONS 


Un  =  1.0  kT  U2i  =  1.5  kT 

Ui2  =  2.0  kT  U22  =  1.5  kT  +  hy 

y  =  500  megahertz 

h  =  Planck's  constant 

»  0.70 
p 12  =  0.30 
Nj  *  1.4  x  101! 

T  =  310°K 

B  =  750  transitions  per  second  per 
milliwatt/cm2  per  molecule 

C2i  =  0.0375  transitions  per  second  per 
molecule 

In  choosing  the  transition  coefficient  B,  a  molecule  with  average 
cross-sectional  area  of  2,500  square  Angstroms  was  considered.  This 
corresponds  to  a  spherical  molecule  with  radius  28.21  Angstroms.  In  a  1 
mW/cm2,  500-MHz  field,  this  molecule  will  be  struck  by  750,000,000  photons 
per  second.  The  coefficient  B  was  chosen  to  provide  one  transition  per 
million  photons. 

Examining  equation  9,  we  see  that  in  the  model  under  discussion  here, 
the  maximum  a  achievable  is  (n21  -  n22)/2.  In  Table  1,  C2i  is  chosen  so 
that  a  perturbation  of  approximately  a/2  is  achieved  at  0.25  mW/cm2. 

Using  equations  4,  5,  and  6,  the  tissue  configuration  prior  to  external 
irradiation  can  be  described.  Specifically,  energy  state  occupancy  and  total 
tissue  energy  can  be  calculated.  These  values  are  shown  in  Table  2  for  the 
model  tissue  at  T  =  310°K. 


P21  =  0.9999 

p22  =  0.0001 

N2  =  1.0  x  106 


TABLE  2.  TISSUE  ENERGY-STATE  OCCUPANCY  AND  TOTAL  ENERGY  WITHOUT  EXTERNAL 

RADIATION 

nu  =  1.2093  x  1011  n2i  =  999900 

n12  =  0.1907  x  1011  n22  =  99.9923 

U  =  6.8069  x  10*10  joules 

(T  =  310°K ) 


In  Figure  2,  tissue  total  energy  is  plotted  as  a  function  of 
temperature.  In  Figure  3,  n22  is  plotted  as  a  function  of  temperature. 
Tables  3  and  4,  respectively,  show  the  data  for  these  figures.  An 
insignificant  increase  in  n22  occurs  as  a  function  of  temperature.  This  is 
due  to  the  closeness  of  the  energy  levels  U2i  and  U22. 


TABLE  3.  TOTAL  TISSUE  ENERGY  AS  A 

FUNCTION  OF  TEMPERATURE 

Temperature  (°K) 

Energy  (joules) 

310 

6.8069  x  10-10 

360 

6.9098  x  10-10 

410 

6.9945  x  10-10 

460 

7.0651  x  10-10 

TABLE  4.  n22  AS  A  FUNCTION 

OF  TEMPERATURE 

Temperature  (°K) 

n22 

310 

99.9923 

360 

410 

460 


99.9933 

99.9941 

99.9948 
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Table  5  shows  the  model  tissue  receiving  0.25  mW/cm2  of  500-MHz 
radiofrequency  radiation.  No  significant  change  in  the  "water"  molecules  is 
noted;  that  is,  no  significant  change  occurs  in  or  nj2*  Neither  is  there 
any  significant  change  in  total  system  energy  U,  indicating  essentially  no 
tissue  temperature  rise.  However,  n22  is  marledly  pumped.  Figure  4  shows 
n2i  and  n22  as  functions  of  irradiation  intensity.  This  figure  indicates 
pumping  at  quite  low  dose  levels;  the  tissue  "enzyme"  is  effectively  very 
"hot"  even  at  0.001  mW/cm2  (recall  Figure  3). 

TABLE  5.  MODEL  TISSUE  UNDER  0.25  mW/cm2  OF  500-MHz  RADIOFREQUENCY  RADIATION 

nn  =  1.2093  x  1011  n21  =  749,959 

ni2  =  0.1907  x  1011  n22  =  250,040 

U  =  6.8069  x  10"!0  joules 
(T  =  310°K) 

A  small  BASIC  language  computer  program  was  written  to  calculate  the 
data  presented  in  Tables  2  through  5  and  in  Figures  2  through  4.  A  listing 
of  this  program  is  provided  in  Appendix  B. 

What  is  the  meaning  of  the  data  presented  in  the  above  discussion? 
Overall,  the  above-outlined  elementary  model  reinforces  the  point: 

"IF  A  MOLECULAR  TRANSITION  IS  SHIELDED  FROM  COLLISIONAL  STIMULATION,  BUT  CAN 
BE  STIMULATED  BY  ELECTROMAGNETIC  RADIATION,  THE  RADIATION  EFFECT  CAN  PUMP  THE 
TRANSITION." 

Two  questions  are  immediate:  (1)  What  may  be  the  biomedical  consequences  of 
pumping  a  transition?  (2)  How  can  we  study  this  phenomenon  experimentally? 

Changing  the  state  of  a  molecule  can  change  its  chemical  reactivity. 
The  biomedical  impact  of  pumping,  then,  can  be  expected  to  shift  biochemical 
reaction  rates  and  pathways  in  the  living  being.  Considering  this  phenomenon 
from  the  biochemical-evolutionary  point  of  view  is  worthwhile.  The  state  or 
conformation  of  a  biomacromolecule  is  very  Important  to  its  reactivity.  As  a 
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Hcjure  4.  Biomacromolecule  energy  level  occupancy  as  a  function  of  exposure 


species  evolves,  key  molecular  conformations  may  be  shielded  from  collisional 
transitioning  to  promote  macromolecular  stability.  These  conformations, 
then,  could  be  vulnerable  to  radiofrequency  radiation. 

At  least  four  research  teams  have  investigated  chemical  reaction  rates 
in  electromagnetic  fields.  Chen  et  al.  (1)  observed  a  small  shift  in  the 
equilibrium  constant  of  an  iron  oxidation-reduction  system  in  a  high  magnetic 
field.  Wayland  and  Brannen  (2)  required  an  electric  field  term  to  explain 
their  thermal  inactivation  curves  for  Bacillus  subtil  is.  Seto  and  Hsieh  (3) 
provided  some  very  interesting  theory  and  experimental  work  with  Escherichia 
col i ,  while  Karube  et  al .  (4)  had  an  intriguing  result  with  membrane-bound 
lipase. 

From  the  medical  point  of  view,  changes  in  chemical  reaction  rate  might 
particularly  be  detected  in  clinical  measures  such  as  tests  of  thyroid  func¬ 
tion  and  of  glucose  metabolism.  Baranski  and  Czerski  (5)  report  research 
along  these  lines. 

The  question  of  pumping  of  biomacromolecules  can  be  examined  experi¬ 
mentally,  or  indirectly,  as  indicated  above,  by  searching  for  biomedical 
evidence  in  animals  or  in  in-vitro  chemical  reaction  experiments.  More 
direct  evidence  can  be  obtained  via  spectroscopic  measurements  and  molecular 
models  as  pursued  by  Prohofsky  and  his  colleagues  {6,  7).  This  research 
group  has  produced  evidence  of  stimulation  of  homopolymer  DNA  and  RNA  mole¬ 
cules  in  electromagnetic  fields  by  gigahertz  radiofrequency  radiation.  The 
biomedical  relevance  of  this  work  needs  further  research.  Of  paticular 
interest  is  the  question  of  whether  molecular  stimulation  can  occur  at  field 
frequencies  of  less  than  1  GHz. 

To  summarize  Part  I  of  this  report,  the  model  tissue  presented  here 
suggests  the  possibility  of  molecular  stimulation,  or  pumping,  by  electro¬ 
magnetic  radiation.  Such  a  phenomenon  could  be  expected  to  lead  to  changed 
chemical  reactivity,  and  research  along  these  lines,  found  in  the  literature, 
is  mentioned.  The  possibility  of  linking  molecular  stimulation  with  biomed¬ 
ical  effects  is  broached.  The  concept  of  molecular  stimulation  or  pumping  by 
electromagnetic  radiation  requires  more  and  deeper  investigation  both  from 
the  point  of  view  of  biophysics  and  from  the  vantage  point  of  biomedical 
effects.  The  present  picture  is  far  from  complete. 


As  a  small  step  toward  an  improved  understanding  of  molecular  stimula¬ 
tion,  a  more  detailed  model  of  macromolecular  responses  in  an  electromagnetic 
field  is  presented  in  the  next  section  of  this  report. 

PART  II:  ON  MACROMOLECULAR  MOTION  IN  AN  ELECTROMAGNETIC  FIELD 

A  very  simple  one-dimensional  chain  molecule  is  used  to  illustrate 
basic  concepts  and  principles.  The  model  molecule  consists  of  n  components 
connected  by  n-1  forces  having  linear  spring  properties.  The  coordinates  for 

the  components  are  xq,  xj,  . .  xn_i;  each  component  has  mass  m,  and  each 

connecting  force  has  force  constant  k.  Each  component  carries  electric 

charge  q0.  qj,  . ,  qn_i*  Figure  5  presents  the  model  schematically.  The 

mathematical  treatment  presented  in  this  section  is  an  adaptation  and  exten¬ 
sion  of  the  analysis  provided  by  Muzushima  (8).  The  model  is  a  simple 
application  of  lattice  theory  (9,  10).  The  setting  is  nonquantum  mechanical 
or  classical  in  approach. 

The  equations  of  motion  for  the  model  molecule  are: 
mx’o  =  -  k(xg  -  xj)  +  qo  E  sin  wt 

mx'i  *  k(xo  -  xj)  -  k ( x i  -  X£)  +  qi  E  sin  u>t 

mx'2  =  k(xj  -  X2>  -  k(x2  -  X3)  +  q2  E  sin  ait  (10) 

•  •  • 

•  •  • 

mx'n-i  =  k(xn_2  -  Xn_i)  +  qn.j  E  sin  wt 

where  E  is  the  peak  electric-field  strength. 

When  divided  by  m  and  rearranged,  these  equations  of  motion  become: 

x’o  =  -lxo  +  lxj  +  eo  sin  wt 

XI  =  Ixq  -  21xi  +  1x2  +  £1  s^n 

x’2  =  Ixi  -  21x2  +  1x3  +  £2  sin  ut  (11) 

•  •  • 

•  •  • 

*n-l  =  1xn-2  -  1  n-1  +  ^n-1  sin  <*>t 
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where  1  =  k/m  and  =  q^  E/m. 

Equations  11  can  be  written  in  matrix  form  as: 


X  = 
— 

-LX  + 

^  sin  uit 

where 

“*0  “ 

“eO  ~ 

X1 

el 

1 

X  = 

x2 

£_  = 

e2 

xn-l 

en-1 

jnd 

1 

-1 

0 

0 

0 

-1 

21 

-1 

0 

0 

0 

-1 

21 

-1 

0 

L  = 

0 

0 

-1 

21 

-1 

Now,  a  matrix  A  Is  sought  such  that 
Ata  =  I  AtlA  =  R 

where  AT  is  the  transpose  of  the  matrix  A,  I  is  the  Identity  matrix,  and 
is  a  diagonal  matrix.  For  the  problem  under  consideration  here,  such  a 
matrix  A  can  be  found.  Specifically,  if 


a00 

a01 

a02 

•  •  • 

a0,  n-1 

a  io 

all 

al2 

•  •  • 

al,  n-1 

a20 

a21 

a22 

•  •  • 

a2,  n-1 

an-l,  0 

• 

an-l,  1 

• 

an-l,  2 

•  •  • 

an-l,  n 

(12) 
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then 

aii  =  (  J-  i1/2  cos  L  ( j+1/2 )  1»/nJ  (13) 

J  n+1 

and  writing  R  as. 


the  rj  are  given  by: 

rf  -  21U  -  cos  (  JJ  (14) 

Equations  13  and  14  are  shown  in  Mizushima  (8). 

Define  =  ri.  The  terms  xi  will  be  very  Important  as  they  represent 
natural  frequencies  associated  with  the  molecular  motion,  as  shall  be  seen 
below. 

Define  new  molecular  coordinates  £,  using  the  matrix  A,  as 
X  =  A£ 

Then  equation  12  becomes 

A£  =  -LA£  +  £  sin  <*>t 

By  multiplying  this  equation  from  the  left  by  AT,  the  following  Important 
equation  is  obtained: 

£  =  -ATLA£  +  At£  sin  “t 
or 

l  »  -RC  +  ATe  sin  wt  (15) 
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Written  In  component  form,  equation  15  becomes  the  set  of  equations 


£o  '  -Vo  *  U0  Vj!  s1"  “* 


£<  -  'V, 


+  l  1  a .  .£ 

j=0  '1 


2  n_1  T 

£  =  -A  £  +  f  )  a 

n-1  n-1  n-1  Ljt0  n-1 

where  aL  Is  the  i,  jth  element  In  AT 


]  sin  wt 


e  1  sin  (^t 

j  j 


Comparing  the  above  set  of  equations  with  the  set  given  before  as 
equation  11,  we  see  that  in  the  above  set  each  equation  involves  just  a 
single  5fj  while  with  equation  11,  each  equation  involves  more  than  one  Xf. 
The  set  of  equations  labelled  equation  16  thus  defines,  in  a  sense,  indepen¬ 
dent  components  of  motion.  In  fact,  the  ^'s  are  called  normal  or  collective 
components,  and  the  modes  of  motion  they  describe  are  called  normal  modes. 

Looking  again  at  equation  16,  we  see  n  independent  oscillators.  The 
solution  can  be  found  by  elementary  methods.  Specifically,  when  xi  * 


£if(t)  =  £f(0)cOS  Aft  +  [CftOj/Af]  sin  Aft 


Ml  el[  “  ]1 


T  u>  sin  U)t  sin  Aft 


j=o  «  *  \2rJ  « 


and,  when  Af  = 


£j(t)  =  CjfOcos  cut  +  [^(O/cuJsin  cot 
n-l  T  1 

+  l  l  a..  e.J[ — _ J [si n  (ot  -  cot  cos  tot]  (18) 

j=0  1J  J  2<o2 

The  values  £^(0)  and  ^.(0)  are  the  coordinate  position  and  velocity, 
respectively,  at  time  t  =  0  which  is  considered  to  be  the  beginning  of  the 
macromolecular  exposure  to  the  electric  field.  These  values  reflect  the 
thermal  state  of  the  molecule  just  prior  to  imposition  of  the  field.  Also, 
note  that  Aj  defines  the  intrinsic  frequency  response  of  the  ith  oscillator. 

The  molecular  perturbation  due  to  the  electric  field  is  symbolized  by 
4|(t),  and  these  perturbations  are  given  by  equation  IP  and  equation  20 
below.  For  Aj  *  to 


r  H- 1  -r 

?  (t)  *  [  l  a  e  ][_!L_ 
i  ijjJLA2.U)2 


]l 


sin  ait 


sin  a^ t 


j=0 


(19) 


and  for  A.j  =  u> 

E  n-l  T  i 

C  ( t )  =  [  l  a  e.][ — —  J [ s i n  u>t  -  wt  cos  u>tj  (20) 

1  j=0  1J  J  zJ 

Note  at  this  point  that  for  Aj  t  w,  the  perturbation  (t)  remains  bounded; 
but  for  a.j  =  a),  s|(t)  grows  linearly  with  time.  It  is  instructive  to  look  at 
the  energy  of  the  molecule  as  a  function  of  the  electromagnetic  perturbation. 
The  kinetic  energy  of  the  molecule  is 

1  »T  •  1  »T  • 

K.E.  =  _mXX=-m5_C 

while  the  potential  energy  of  the  molecule  is 

1  T  1  T 
P.E.  =  _mXLX=-m^Rj, 
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Using  the  above  facts,  we  state  the  molecular  energy  due  to  electromagnetic 
perturbation  of  the  ^  coordinate  as 


•2122 
S.  +  -  mx.  £. 

i  2  1  7 


Using  this  last  equation  and  equations  19  and  20,  we  find  that  for  Xi  *  u> 

1  2222  n*1T  2  2 

W=  —  mlai  /( x  -w  )  II  )  a .  .e  1  jl+cosajt  + 
i  2  i  j_Q  ij  j 

9  7  7 

{ X(/ur)sinewt  -  2  cos  wt  cos  x^t  -  2(x^/ui)sin  wt  sin  x^tj  (21) 
and  for  X,-  =  u 

2  n“*  T  2  2  2 

W.  =  (m/8ui  )[  l  a..e.J  11  +  ut  |  (22) 

1  j=0  1J’  j 

These  last  two  equations  are  of  particular  interest.  When  the  stimulating 
frequency  m  is  not  equal  to  a  natural  frequency  of  the  molecule  Xj ,  the 
absorbed  energy  is  a  finite  quantity,  as  indicated  by  equation  21.  This 
absorbed  energy  can  be  compared  with  the  thermal  energy  in  an  individual 
normal  mode,  which  is  kT  (k  =  Boltzmann  constant,  T  =  temperature  in  °K)  by 
the  equivalence  principle  (11).  When  the  stimulating  frequency  w  equals  a 
natural  mode  frequency  of  the  molecule  (the  resonance  case),  the  molecule 
gains  energy  quadradical ly  with  time,  as  seen  in  equation  22.  If,  at  reso¬ 
nance,  the  molecule  cannot  adequately  dissipate  absorbed  energy,  molecular 
disruption  would  be  expected. 

Based  upon  equations  21  and  22,  but  particularly  equation  22  with  the 
resonance  response,  developing  a  fuller  understanding  of  molecular  dynamics 
in  a  radiofrequency  field  is  clearly  important.  Two  questions  are  immediate. 
(1)  What  are  the  natural  frequencies  of  biomacromolecules  in  tissue?  (2) 

What  energy  dissipative  mechanisms  are  available  to  biomacromolecules  in 
tissue?  These  two  questions  will  now  be  addressed  in  the  context  of  the 
simple  linear  chain  model  already  under  discussion. 
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B I OMOLEC ULAR  NATURAL  VIBRATIONAL  FREQUENCIES 


Regarding  the  natural  frequencies  of  biomacromolecules  in  tissue,  equa¬ 
tion  14  can  be  used  to  provide  a  first  approximation.  Using  k  =  3  millidynes 
per  Angstrom,  m  =  200  gram  molecular  weight,  and  N  =  2000,  the  lowest  ten 
natural  frequencies  of  the  chain  model  are  given  in  Table  6.  The  value  of  k 
(3  millidynes  per  Angstrom)  reflects  the  type  of  force  constant  one  can  find 
in  a  carbon-carbon  bond.  The  m  and  N  values  employed  in  constructing  Table  6 
provide  a  macromolecule  with  gram  molecular  weight  equal  to  400,000. 

TABLE  6.  THE  LOWEST  TEN  NATURAL  FREQUENCIES  IN  A  MODEL  MACROMOLECULE 

GHz 

0.000 

7.512 

15.025 

22.537 

30.050 

37.562 

45.075 

52.587 

60.100 

67.612 

In  Table  6  the  lowest  frequency,  zero,  corresponds  to  direct  accelera¬ 
tion  or  translation  of  the  entire  molecule  by  the  electromagnetic  field.  The 
next  lowest  frequency  is  an  internal  vibratory  mode.  This  mode  and  several 
others  are  in  the  low  gigahertz  range.  If  the  model  were  made  longer  and/or 
heavier,  natural  vibratory  modes  could  be  found  in  the  100  to  900  MHz  range. 
Such  molecules  would  be  quite  large  linear  structures,  however.  More 
complete  calculations  and  spectroscopic  data  are  needed  to  determine  the 
extent  to  which  vibratory  resonances  occur  below  1  GHz.  However,  a  brief 
glance  at  most  spectroscopy  texts  will  rapidly  convince  the  reader  that 
rotational  modes  are  commonly  found  below  1  GHz.  Again,  future  research  will 
determine  the  occurrence  and  biomedical  meaning  of  these  modes. 
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AVAILABLE  DISSIPATIVE  MECHANISMS 

Two  d1ss4vjtive  mechanisms  come  immediately  to  mind  for  consideration: 
radiation  of  energy,  and  loss  of  energy  via  collisional  events. 

The  molecular  charges  q,,  when  they  move,  will  lose  energy  by  setting 

up  a  radiation  field.  A  damping  force,  Fj,  is  associated  with  this  radiation 

as  the  moving  charge  interacts  with  its  own  field.  This  force  is  given  by 

,2  2  i  ..  2  3  ...  ,  , 

F  =  (q  /8ue  c  r  )  x  -fq  /6tte  c  )  x  +  •••♦  (23) 

i  1  i  0  O'  i  l  0  i 

In  this  equation,  c  is  the  speed  of  light,  is  the  capacitivity  of  the 
vacuum,  and  rg  is  the  radius  of  the  ith-charged  particle.  Higher  order  terms 
in  this  force  equation  are  neglected  for  they  converge  towards  zero  as  rg 
approaches  zero  (8).  Equation  23  will  be  incorporated  into  equation  10. 

Collisional  events  will  exert  a  damping  force  on  molecular  motion  also. 
In  general,  the  collisional  damping  on  the  ith  molecular  component  will  be 
3^  where  3^  is  a  coefficient  of  friction. 

With  the  above  considerations,  equation  10  becomes: 

••  2  ••  2 . . . 
mx  =  -k ( x  -x  )  -  a  q  x  +  bq  x  -  3  x  +  q  E  sin  u>t 
0  0  1  0  0  0  0  0  0  0  0 

••  2  ••  2  .  .  . 

mx  =  k(x  -x  )  -  k(x  -x)-aq  x  +bq  x  -  3  x  ^  q  E  sin  u>t 

1  0  1  12  11111111 

••  2  ••  2. . . 

mx  =  k(x  -x  )  -  k ( x  -x  )  -  a  q  x  +  bq’  x  -  3  x  +  q  E  sin  at  (24) 

2  12  23  2  2  2  2  2  22  2 


2  ..  ? 

mx  =  k ( x  -x  )-a  q  x  +  bq"  k‘  -3  x  +  q  E  sin  wt 
n-1  n-2  n-1  n-ln-1  n-1  n-1  n-1  n-1  n-1  n-1 
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where  a^.  =  and  b  =  Rearranging,  we  obtain  the 

following: 

-bVo  *  l""aoqo)xo  *  eoxo  '  "kx0  4  kxi  4  "oE  s1n  “t 
2. . .  2 

-bq  'x'  +  (m+a  q  jx  +  v  x  *  kx  -  2kx  +  kx  +  q  E  sin  cot 

11  11111  0  1  2  1 

2 . . .  2  •• 

-bq2'x2  +  (m+a2q2Jx2  +  6^  =  kx^  2kx?  +  kx3  +  q^E  sin  cot  (25) 


2  ... 


■b1n  /’  ,  +  q  Jx  -  +  B  ,x  =  kx  -  kx  ,  +  q  E  sin  cot 

n-1  n-1  n-1  n-1  n-1  n-1  n-1  n-2  n-1  n-1 


Let  Tq  ro  -  r0  “  *  •  •  =  r0  ^  =  r0’  so  a0  =  al  =  a2  =  •••  =  an-l  =  a» 
and  let  qjj  =  qf  =  q|  =  ...  =  q^  =  q* ,  6Q  =  6j  =  62  =  ...  =  6^  =  e,  in  = 
m  +  aq2,  x  =  -bq2/  m,  1  =  k/m,  6  =  6/m,  and  =  q^/m.  Then  equation  25 
becomes: 


t  x  +  x  +  6  x  =  -lx  +  lx 
0  0  0  0  1 


+  e  sin  cot 
0 


x  x  +  x  +  6  x  =  lx  -  21  x  +  lx  +  e  sin  cot 
11  1  0  12  1 


x  x?  +  x2  +  6  x^  =  lx:  ‘  21  x2  +  lx3  +  e2  sin  cot 


(26) 


T  xn  1  +  Xn  1  +  yx  ,  =  1x  o  -  1x  ,  +  e  sin  cot 

n-1  n-1  n-1  n-2  n-1  n-1 


In  matrix  form  equation  26  is: 


xX  +  X  +  6  X  =  -l  X  +  e  sin  cot 
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Again  using  X  =  H,  we  obtain 

t  i +  £  +  6_c  =  -Rc  +  AT£  sin  tot 

where  R  is  a  diagonal  matrix  with  entries 

.  .  in  .2 

r.  =  2f 1 1  -  cos  ( — ]J  =  x. 
i  1  vn  J  l 

The  ith  component  of  the  matrix  equation,  equation  27,  is 


2  n-1  T 

t  £  +  £  +  +  x.  5.  =  >  a. .el  sin  wt 

11  111  j^Q  1J  j" 


(27) 


(28) 


Solving,  we  find  that  for  x,-  *  to,  the  molecular  perturbation  due  to  the 
electric  field  of  the  ith  coordinate  is: 


S 


E 

i 


(t) 


n-1 

L 

j=0 


.22.  .  .  3^ 

(Xj-to  J  Sin  tot  -  (gio-Tto  J  COS  <ot 

.2  2  2  7  r~3~2 

(X^-to  J  +  (pio-Tto  J 


When  x-j  =  to  ( i . e . ,  at  resonance) 


cos  <ot 

(gtO-TtO^) 


n"X  T  . 

i  1  a . -e . 

l  lj=0  u  J 


and 


L  = 


"-1  T  . 

1  3 . .G. 

j=0  ’J  J 


to  si  n  tot 

(010-TtO^) 


Using 


2  2  2 
W  =  (l/2)m£  +  (l/2)mx  £ 
i  i  ii 


we  find  that 


W  =  (1/2)  m 

i 


) 


2 


(29) 


(30) 


(31) 


(32) 
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From  this  last  equation  we  see  that  molecular  energy  remains  finite  at 
resonance  when  we  consider  the  effect  of  dissipative  mechanisms  (remember 
that  t  is  a  negative  number). 

With  equations  28  and  32  in  hand,  we  can  appropriately  reconsider  the 
400,000  molecular  weight  macromolecule  studied  earlier.  Recall  that  this 
macromolecule  was  defined  by  the  following  parameters:  k  =  3  millidynes  per 
Angstrom,  m  =  200  gram  molecular  weight,  N  =  2000.  For  the  present  analysis, 
we  will  consider  charges  symmetrically  placed  on  the  molecule.  Specifically, 
let  us  assume  that  qg  through  qg65  equal  0.0015015  times  the  charge  on  one 
electron  (-1.6021917  x  10_19  Coulomb),  qgsg  through  qi333  equal  0.0015015 
times  the  charge  on  one  positron,  and  qj334  through  q^gg  equal  0.0015015 
times  the  charge  on  one  electron.  This  charge  distribution  provides  one 
electron  charge  spread  over  the  left  and  right  outer  one-third  parts  of  the 
chain  model  molecule,  and  provides  one  positron  charge  over  the  central 
one-third  portion  of  the  molecule. 

For  rg,  we  use  value  2.817939  x  10-15  m,  a  calculated  electron  radius 
(8).  Of  course,  eg  =  8.85419  x  10"12  Farads/m,  c  =  2.9979  x  108  m/second. 

The  only  term  left  undecided  is  the  friction  coefficient  g.  In  Appendix  C  we 
see  that 

B  =  "72t  (33) 

where  m  is  the  mass  of  the  macromolecular  component  and  t  is  the  time  con¬ 
stant  of  the  molecule  relative  to  energy  dissipation  by  colli sional  mecha¬ 
nisms. 

Working  first  with  equation  28,  we  see  immediately  that,  since  m  is 
only  very  slightly  larger  than  m,  the  molecular  natural  frequencies  with 
dissipation,  Xi»  are  essentially  unchanged  from  the  natural  frequencies  Af 
without  damping.  Working  at  the  lowest  frequency  of  7.512  GHz  and  using 
equation  32,  we  calculate  the  data  in  Figure  6.  Appendix  D  provides  the 
computer  program  used  to  accomplish  these  calculations. 

Studying  Figure  6,  we  again  see  a  phenomenon  that  was  first  noted  in 
Part  I  of  this  report.  Namely,  when  the  molecular  mode  is  shielded  from 
collisional  exchanges,  the  electromagnetic  field  is  able  to  significantly 
excite  the  mode.  Note  specifically  that  at  higher  time  constants,  mode 
energies  considerably  above  kT  are  seen.  Further  analysis  and  experimenta¬ 
tion  are  needed  to  determine  whether  any  such  high  time  constants  actually 
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Figure  6.  Normal  mode  energy  vs.  RFR  exposure  and  molecule  relaxati 


obtain  in  living  tissue.  As  discussed  in  Part  I,  molecular  stimulation  or 
pumping  could  be  detected  as  changed  chemical  reactivity  or  by  spectroscopic 
measurements. 

Figure  7  presents  a  schematic  of  molecular  structure  which  might 
evidence  modes  of  motion  that  are  shielded  from  colli  si onal  exchanges  but  are 
susceptible  to  electromagnetic-field  stimulation.  The  inner  molecular 
moieties  are  proposed  as  the  potentially  susceptible  moieties.  In  Figure  8 
the  molecular  structure  introduced  in  Figure  7  is  shown  interacting  with 
another  compound  to  form  a  more  complex  entity.  This  complex  entity  is 
intended  to  represent  an  antigen-antibody  combination  or  another  similar 
receptor-site  interaction.  Again,  the  inner  molecular  moieties  would  be 
shielded  from  collisional  energy  exchanges  but  could  be  potentially  suscepti¬ 
ble  to  perturbation  by  electromagnetic  fields.  This  concept  that  receptor- 
site  mechanisms  could  be  perturbed  by  electromagnetic  radiation  leads  to  the 
conjecture  that  biological  events  that  involve  coding  and  decoding  may  be 
susceptible  to  radiofrequency  radiation.  For  example,  DNA-RNA  transcription 
events  could  be  vulnerable,  as  could  be  antigen-antibody  immunological 
responses.  Autoimmune  antibody  has  been  noted  in  animals  chronically  irradi¬ 
ated  with  radiofrequency  radiation.  This  leads  us  to  consider  whether  radio¬ 
frequency  radiation  might  not  be  implicated  in  the  genesis  of  collagen 
diseases  such  as  lupus  erythematosus  or  scleroderma  (12). 

PART  III:  DISCUSSION 

In  Part  I  of  this  report  an  elementary  statistical  physics  model  of 
tissue  was  presented.  This  model  suggested  that  molecular  transitions  which 
were  shielded  from  collisional  energy  exchanges  could  be  stimulated  by 
electromagnetic  radiation.  Part  II  of  this  report  presented  a  simple  model 
linear  chain  macromolecule.  This  macromolecule  exhibited  natural  frequencies 
in  the  low  gigahertz  range.  Once  again  it  was  seen  that  if  a  molecular  mode 
were  shielded  from  collisional  energy  exchanges,  electromagnetic  stimulation 
could  obtain.  Biomedical  consequences  of  macromolecular  electromagnetic 
stimulation  would  be  effects  following  changed  macromolecular  chemical 
reactivity,  as  stimulation  of  molecular  modes  of  motion  could  potentially 
change  the  ability  of  the  molecule  to  engage  in  specific  chemical  reactions. 
The  pivotal  point  now  appears  to  be  the  question  of  how  a  macromolecule 
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dissipates  energy  absorbed  from  an  electromagnetic  field.  A  central  question 
concerns  the  collisional  mechanisms  available  to  a  macromolecule  in  living 
tissue.  This  question  can  be  examined  by  studying  chemical  reactions  in 
tissue  or  tissue-equivalent  material.  Also,  more  directly,  electromagnetic 
absorption  data  obtained  from  macromolecules  in  a  tissue  milieu  will  help 
elucidate  dissipation  mechanisms. 

In  this  report's  treatment  of  a  biomacromolecule  imbedded  in  a  medium, 
electromagnetic  radiation  that  will  be  incident  upon  the  molecule  from  adja¬ 
cent  macromolecules  as  well  as  from  external  sources  is  not  included.  This 
Internal  radiation  exchange  needs  and  warrants  further  investigation. 
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APPENDIX  A 


DERIVATION  OF  EQUATIONS  4  AND  5 

Molecules  of  two  types,  #1  and  #2,  are  treated.  Type  #1  molecules  can 
occupy  energy  levels  U^,  U12,  U13,  ....  U1L,  and  type  #2  molecules  can 
occupy  energy  levels  U2j,  U22,  U23,  ....  U2M.  The  symbol  x^j  represents  the 
number  of  molecules  of  type  i  in  energy  level  j.  To  determine  the  most 
likely  values  n^,  of  the  x^'s,  the  probability  of  the  different  possible 
system  configurations  is  maximized  under  the  constraints  of  (a)  constant 
total  molecules  of  each  type,  and  (b)  constant  total  system  energy.  Thus, 
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is  maximized  with  constraints 


L 

l  xij  =  Ni  =  total  number  of  type  #1  molecules 
j=l 

M 

l  x2j  =  N2  =  total  number  of  type  #2  molecules 
j=l 

L  M 

l  x^Ujj  +  l  XojUo..  =  U  =  total  system  (mixture)  energy 
j*l  J  J  j=l  J  J 

The  function  P(x^,  ...;  x^L,  x2p  ...,  x^)  is  the  product  of  two 

multinomial  distributions  where  the  constants  p^j  are  the  probabilities  that 
a  molecule  of  type  i  can  enter  its  jth  energy  level  (1). 
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The  maximization  is  achieved  by  working  with  the  natural  logarithm  of  P 
using  Lagrange  multipliers  a,  b,  and  c,  and  solving  the  following  equation 
(2): 

dlnP  +  adNj  +  bdN2  +  cdU  s  0 
Note  that 

dS  =  kdlnP 

where  S  is  system  entropy  and  k  is  the  Boltzmann  constant;  and  recall  the 
thermodynamic  equation 

TdS  +  ujdNj  +  u2dN2  -  dU  =  0 

to  determine  immediately  that  (3) 

c  =  -(kT)'1 

Using  Stirling's  formula  (In  R!  «  RlnR  -  R), 

L  M 

dlnP  *  l  (In  pXj  -  In  x^)  dx^  +  l  (In  - 

j=l  j=l 

whence 

dlnP  +  adN^  +  bdN2  +  cdU  «• 

L 

l  (Inpjj  -  lnxjj  +  a  +  cUjjldxjj 
j=l 


In  Xytdny 


<A2 


+  l  Hnp2j  -  lnx2j  +  b  +  cU2j)dx2j 

j=l 

and,  for  occupancy  numbers  x^j,  such  that  Stirling's  formula  is  accurate, 

P ( xi l ,  xj2,  x2i,  ...,  x2m)  is  maximized  by  xij=nij  and  x2j=n2j»  where 


nlj  =  plj 


e+ae-ulj/kT 


and 


"2J  =  P2je*VU2"kT 


34 


Using  the  relations 


\ «1J  ■  N1 
J=1 

M 

l  x2j  =  N2 
J=1 

the  constants  a  and  b  can  be  found,  providing  equations  4  and  5  of  the  text. 
Differentiating  equation  A2  confirms  that  the  n^j's  maximize  P. 
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APPENDIX  B 


BASIC  LANGUAGE  COMPUTER  PROGRAM  FOR  CALCULATING  THE  DATA 
FOUND  IN  TABLES  2  THROUGH  5  AND  IN  FIGURES  2  THROUGH  4 

10  REM  THE  TITLE  OF  THIS  PROGRAM  IS  EQUIL 
20  DIM  Pl(2) ,  P2(2) ,  Ul(2) ,  U2(2) 

30  DIM  A1 { 2 ) ,  A2{2) 

40  P 1 ( 1 ) =0  -  70 

50  Pl(2)=0.30 

60  P2(l)=0.9999 

70  P2 ( 2 ) =0 . 0001 

80  Nl=l -4E11 

90  N2=1.0E6 

100  T=310 

110  V=500E6 

120  H=6.6252E-34 

130  K=1 .3804E-23 

140  Ul(l )=1 .0*K*310 

150  U1(2)S2.0*K*310 

160  U2(l)=l .5*K*310 

170  U2(2)=1.5*K*310  +  H*V 

180  Z1=P1( 1)*EXP( -U1(1)/{K*T) )  +  P1(2)*EXP(-U1(2)/(K*T) ) 

190  Z2=P2(1 )*EXP(-U2( 1 )/(K*T) )  +  P2(2)*EXP(-U2(2)/(K*T) ) 

200  FOR  I=1T0  2 

210  Al( I )=P1( I )*N1*EXP( -Ul( I )/(K*T) )/Zl 
220  A2( I )=P2( I)*N2*EXP(-U2( I )/{K*T) )/Z2 
230  NEXT  I 

240  E=A1 ( 1 ) *U 1 { 1 ) +A1 ( 2 ) *U 1(2) +A2 { 1 ) *U2 ( 1 ) +A2 ( 2 ) *U2 ( 2 ) 

250  PRINT 

260  PRINT  "TEMP",  T 
270  PRINT 

280  PRINT  ”N1(1)\A1U) 

29f  .RINT  "N1(2)",A1{2) 

300  PRINT  "N2{1)",A2(1) 
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310  PRINT  "N2(2)",A2(?' 

320  PRINT 

330  PRINT  "ENERGY",  E 
340  DIM  C2(2) 

350  8*750 

360  C2(l)=0.0375 

370  C2(2)=C2( 1)*(A2( 1 )/A2(2) ) 

380  PRINT 

390  PRINT  "C2(1)",C2(1) 

400  PRINT  "C2(2)",C2(2) 

410  I=.001 

420  D=(B*I*(A2( 1)-A2(2) ) )/(2*B*I+C2(l)+C2(2) ) 

430  PRINT 
440  PRINT  "I", I 
450  PRINT  "DELTA", D 
460  A2(l )=A2( 1 )-D 
470  A2(2)=A2(2)+D 
480  PRINT 

490  PRINT  "NEW  N2( 1 ) " ,A2( 1 ) 

500  PRINT  "NEW  N2(2)",A2(2) 

510  E=A1(1)*U1(1)+A1(2)*U1(2)+A2(1)*U2(1)+A2(2)*U2(2) 
520  PRINT 

530  PRINT  "NEW  E",E 
540  END 
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DERIVATION  OF  EQUATION  33 


It  is  well  known  that  drag  forces  on  a  surface  are  proportional  to  the 
velocity  gradient  perpendicular  to  the  surface  (1).  Thus  a  term  is 
inserted  in  equation  24.  If  is  the  drag  force,  then  must  be  the 
rate  at  which  energy  is  lost  via  these  forces.  That  is: 


dE  = 
dt 


(Cl) 


Drag  forces  are  the  result  of  collisional  events  between  the  macromolecule 
and  its  solvent.  In  colliding  with  the  macromolecule,  solvent  molecules  can 
change  the  partitioning  of  energy  in  the  macromolecule.  Specifically, 
collisions  can  change  the  energy  of  vibration  or  the  energy  associated  with 
electron  motion  in  the  macromolecule.  Thus  E  may  be  written  as: 

E  “  Ey  +  E0 

where  Ev  is  the  energy  of  vibrational  motion  and  Ee  is  the  electron  energy. 
Equation  Cl  can  be  written 

dE  =  E  =  _Ekv 
dt  m  kv  T 

where  Ej<v  is  the  instantaneous  kinetic  energy  of  vibration  and  T  is  the 
desired  time  constant. 
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APPENDIX  D 


BASIC  LANGUAGE  COMPUTER  PROGRAM  USED  TO 
CALCULATE  DATA  IN  FIGURE  6 


10  REM  CALCULATION  OF  W 
20  SELECT  R 

30  REM  CALCULATE  SUM  OVER  J  OF  A(J I )*Q( J ) 
40  1=1 
50  N=2000 
60  S1=0.0 

70  Q=1.6021917E-19 

80  FOR  J=0  TO  665 

90  F0=-0. 0015015 

100  Z=COS( ( I+0.5)*J*#PI/N)*F0*Q 

110  S1=S1+Z 

120  NEXT  J 

130  FOR  J=666  TO  1333 

140  F0=0. 0015015 

150  Z=COS( ( I+0.5)*J*#PI/N)*F0*Q 

160  S1=S1+Z 

170  NEXT  J 

180  FOR  J=1334  TO  1999 

190  F0=-0. 0015015 

200  Z=C0S((I+0.5)*J*#PI/N)*F0*Q 

210  S1=S1+Z 

220  NEXT  J 

230  PRINT  "SUM(A(JI)*Q(J))\S1 

240  M=200 

250  M*M/6.02E26 

260  REM  B  IS  FRICTION  COEFFICIENT 
270  INPUT  "T",T 
280  B=0.5*M/T 
290  E0=8.85419E-12 
300  C=2.9979E8 
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310  R0*2 . 81 7939E- 1 5 
320  A0*8*#P  I  *E0*C',2*R0 
330  A0*1/A0 
340  B0*6*#P I *E0*C~3 
350  B0=1/B0 
360  E=61.40 
370  L=2*#PI*7.512E6 
380  Ml=M+A0*(F0*Q)/^2 
390  B1=B/M1 

400  T1=-(B0*(F0*Q)^2)/M1 
410  0*1/(B1-T1*L^2) 

420  S2=SQR(2/(N-1) )*$1 
430  S2»(E/M1)*S2 
440  W=0.5*M*(D*S2p2 
450  PRINT  "T",T 
460  PRINT  "W\W 
470  END 


40 


